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Abstract 

The long range character of the exchange coupling between localized magnetic moments indirectly 
mediated by the conduction electrons of metallic hosts can play a significant role in determining 
the magnetic order of low-dimensional structures. Here we consider how this indirect coupling 
influences the magnetic alignment of adatoms attached to the walls of carbon nanotubes. A 
general expression for the indirect coupling in terms of single-particle Green functions is presented. 
Contrary to the general property that magnetic moments embedded in a metal display Friedel-like 
oscillations in their magnetic response, calculated values for the coupling across metallic zigzag 
nanotubes show monotonic behaviour as a function of the adatom separation. Rather than an 
intrinsic property, the monotonicity is shown to reflect a commensurability effect in which the 
coupling oscillates with periods that coincide with the lattice parameter of the nanotube host. 
Such a commensurability effect does not dominate the coupling across semiconducting zigzag or 
metallic armchair nanotubes. We argue that such a long-range character in the magnetic interaction 
can be used in future spintronic devices. 
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I. INTRODUCTION 



More than a decade after the discovery of carbon nanotubes, these nanoscale cylindrical 
structures are still the subject of intensive scientific research due to their intriguing physical 
properties. Significant progress has been made to explain the intrinsic properties of nan- 
otubes but in order to expand the applicability of those systems we need to understand how 
they are affected by the interaction with other objects. Nanotubes interacting with magnetic 
foreign objects are now in focus due to the possibility of implementing the technologically 
promising area of spintronics in molecular structures. In fact, the ability to produce size- 
able changes in the conductance of a nanotube due to an applied magnetic field is one of 
the driving forces in the research of magnetic properties of carbon-based structures.— The 
existence of such a spin-valve effect has potential applications in the development of faster, 
smaller and more efficient nanoscopic magneto-electronic devices. 

Due to their inherent spin asymmetry, the interaction with magnetic foreign objects 
is likely to cause a spin-dependent response on the transport properties of the combined 
structure. Substrates^, substitutional impurities 4 , adsorbed atoms^ and nanoparticles 7 -, 
are some of the different magnetic foreign objects that can interact with carbon nanotubes. 
Among those, transition-metal magnetic adatoms have been reported to produce noticeable 
changes in the spin-dependent electronic structure of carbon nanotubes^ 5 ' 7 

Besides establishing how magnetic impurities affect the electronic structure of the nan- 
otube, it is crucial to understand the nature of the coupling between nearby adatoms since 
the transport properties of a magnetically-doped structure depend on how the impurity mo- 
ments are oriented. Dipolar and exchange interactions are the basic mechanisms defining 
the alignment of the adatoms.— The former decays rather quickly as the moments are moved 
apart whereas the latter depends on both the dimensionality and on the nature of the inter- 
action. Direct exchange coupling results from the overlap between wave functions centred 
at the magnetic impurities but also decays abruptly as the distance between the impurities 
increases. Of indirect nature, the exchange interaction between magnetic impurities medi- 
ated by the conduction electrons of the non-magnetic host is known to decay more slowly 
and plays an important role in determining the overall magnetic alignment of the system. 

The coupling between magnetic moments mediated by the conduction electrons of non- 
magnetic materials has been widely studied, both theoretically and experimentally. The 
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main focus of recent work has been on metallic multilayers^, where magnetic slabs are sep- 
arated by layers of non-magnetic metals. These systems are typically tri-dimensional, with 
translational symmetry preserved in two of the three dimensions. Many of the features of 
the exchange coupling in these systems are directly linked to their geometrical properties 
and therefore are expected to change appreciably for systems of lower dimensionality such 
as carbon nanotubes. In particular, it has been predicted that the indirect exchange cou- 
pling through a one-dimensional electron gas is long-ranged. 10 If this prediction holds for 
the coupling of magnetic adatoms mediated by nanotubes, this long-ranged magnetic inter- 
action could lead to a degree of correlation in the electronic potential of magnetically doped 
nanotubes. From the perspective of fundamental science, such a spin-dependent correlated 
disorder is certain to affect the magneto-transport response of these structures, justifying 
the interest of our study. On the applied side on the other hand, our study is justified by 
the fact that detailed knowledge and, perhaps, control of this magnetic coupling opens the 
possibility of building quantum logic gates in solid state environments, reminding us of the 
huge potential of these structures as future components of spintronic devices. 

With this motivation, in this article we investigate the indirect magnetic coupling between 
two magnetic adatoms attached to the walls of a carbon nanotube. By studying the energetic 
balance of the system as a function of the relative angle between the atomic magnetizations 
we can estimate the nature and magnitude of the coupling across the nanotube. In this way, 
we can obtain not only how the magnetic moments are aligned but also how it decays as 
the adatoms are moved further apart. The sequence adopted in this article is the following. 
We start by deriving a closed-form expression describing the energy required to rotate the 
relative magnetizations of the two adatoms. Since this expression involves matrix elements 
of the single-particle Green function of the system, it is worth presenting in the subsequent 
section a general way of evaluating these matrix elements for achiral carbon nanotubes. 
Results and discussions are presented in section IIV[ followed by conclusions. 

II. INDIRECT MAGNETIC COUPLING 

We consider two magnetic atoms, labelled A and B, adsorbed onto the walls of an in- 
finitely long carbon nanotube and schematically represented in Figure ^ Magnetism in these 
atoms is driven by an intra-atomic Coulomb interaction that, when treated in mean- field 
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FIG. 1: (color online) Schematic representation of the magnetic adatoms on a carbon nanotube. 

approximation through a self-consistent procedure, can be described by an effective spin- 
dependent potential located at the atomic positions. In this way, the electronic structure of 
the entire system is well described by a single-particle Hamiltonian in a basis of localized 
atomic orbitals. In such a basis, the tight-binding-like Hamiltonian is fully determined by 
the on-site potentials and hopping integrals. 

We start by assuming that the magnetic moments of the individual adatoms are initially 
parallel, hereafter referred to as the ferromagnetic (FM) configuration. In this configuration, 
the Hamiltonian of the entire system written in the basis \j) of atomic orbitals centred at a 
site j is given by H = H nt + Ha + H b + V c , where H NT = 12jj> is the Hamiltonian of 

the individual nanotube, Ha = \a)e a (a\ and Hb = \b)eb(b\ are the Hamiltonians associated 
with the isolated atoms A and B, respectively, and Vc = h{\a)(0\ + \0)(a\ + \b)(n\ + \n){b\} 
refers to the coupling between the adatoms and the nanotube. The parameters 7, e a , 
and t are all matrices in spin and orbital spaces and correspond to the hopping between 
nearest-neighbour sites in the nanotube, the on-site potentials of atom A, of atom B and the 
hopping between the nanotube atoms and the adatoms, respectively. Likewise, the basis \ j) 
represents vectors in the same linear space. It is evident from the expressions above that 



sites j = a and j = b label the two adatoms and j = and j = n label the nanotube sites 
that are coupled to the magnetic atoms. 

Since we are interested in evaluating the energy required to rotate the magnetic moment 
of one adatom relatively to the other, it is instructive to start by expressing the total energy 
of the system in the FM configuration, i.e., 



S 



FM 



I _|_ e 0(ui-n) 



ImTr^G,^ 



(1) 



where j3 = l/ksT, kg is the Boltzmann constant, T is the temperature, /i is the Fermi level, 
Gjj(u) is the Green function of an electron with energy u localized at site j and Tr stands 
for the trace over the orbital and spin degrees of freedom. To rotate the magnetic moment 
of one adatom, say atom A, we introduce a perturbation V(9) localized on this atom that 
describes how the effective potential is altered by the rotation of an angle 9. In terms of the 
Pauli matrices a x , a y and a z , it is given by V{9) = —V x ® [(cos9 — 1) o z + sin 9a x ], where V x 
is a matrix in orbital space representing the strength of the local exchange potentials. The 
energy change that results from the perturbation V(8) can be expressed by the respective 
changes produced in the Green function Gjj summed over all sites j (including j = a and 
j = b). Making use of suitable sum rules for Green functions, we can write this energy 
change a a 11 ' 12 



AS (6) 



du 



ImTrln[l-G o , a (u;)F(0)] 
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where G a ^ a represents the Green function in the FM configuration projected at the atomic 
site j = a. Further manipulation of Eq. (j2J) leads to an explicit dependence of the integrand 
on the spin indices, i.e., 
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(3) 
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where G^^uj) represents the propagator between sites j = £ and j = m for electrons of 
spin a and energy u, and the trace is now only over orbital indices. Since the cos 9 in the 
integrand above is the only ^-dependent term, it is reasonable to assume a Heisenberg-like 
angular dependence for AS and conclude that the minimum energy is either at 9 = or 
9 = 7i. By defining the coupling as AS(9 = tv) we can infer not only about the strength 
of the indirect interaction between the adatoms but also whether their moments display a 
parallel or anti-parallel alignment. 



For the sake of simplicity the electronic structure of the system will be here treated 
within the single-band tight-binding model. The expressions above are very general and by 
no means restricted to such a simple case. The results here obtained can be easily extended 
to a multi-orbital description but bring no qualitative difference. Having transition-metal 
atoms in mind, the adatoms are described by a 5-fold degenerate <i-band with the appropriate 
occupation to represent typical magnetic materials. The carbon nanotube is also known to be 
well described by a single-band tight-binding model with suitable parameters that sucessfully 
reproduce the experimental data for the nanotube electronic structure. In this way, rather 
than matrices in orbital indices, all quantities in the integrand of Eq. (JHJ) become scalar. 

It is worth mentioning that the dependence of AS on the distance between the two 
magnetic adatoms is contained entirely in the matrix elements G a ah and G% a . To make this 
dependence more explicit we can express the Green function G associated with the system 
in the FM configuration in terms of another Green function associated with the separate 
parts (also in the FM configuration). In other words, by considering the coupling h between 
the nanotube and the adatoms as a perturbation, we can use Dyson's equation to write the 
matrix elements of G in terms of those associated with their isolated components. It is then 
straightforward to see that 



where C = , g ^ ak '-J^ , P = Gbbh 2 G nn and Q = , r h ' h ui r — • Since the Green function Q 



represents the electronic propagator associated with the nanotube and adatoms in isolation, 

it is evident that Q aa = and Qbb = ■ Regarding the Green functions Q 0n and Q n o, 

these are the only matrix elements associated with impurity-free nanotubes that carry the 
spatial dependence of the indirect magnetic coupling between the adatoms. The indirect 
character of the coupling is evident in the expressions above when written in terms of the 
electronic propagators between the sites that are coupled to the magnetic adatoms. Put in 
another way, the presence of those off-diagonal propagators in the expression for the coupling 
indicates that the magnetic information between the moments is being carried back and forth 
by the conduction electrons of the nanotube host. Since these propagators play such an 
important role in determining the most energetically favourable magnetic configuration, in 
the following section we present a fully analytical method providing closed-form expressions 
for arbitrary matrix elements of the nanotube Green function, among which the elements 
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Q 0: n and Q nfi . 



III. GREEN FUNCTION OF CARBON NANOTUBES 

Carbon nanotubes are fullerene-like structures that can be regarded as graphene sheets 
wrapped up in cylindrical shape. The electronic structures of both graphene and carbon 
nanotubes are well described by a tight-binding model for the 7r-orbital, the only difference 
between them being the wave vector quantization along the circumferential direction of the 
tubes. With such a simple dispersion relation describing the band structure of pure tubes, 
one can obtain the corresponding single-particle Green functions without having to evaluate 
them numerically. In this section we derive an analytical closed-form expression for a general 
matrix element of the single-particle Green function of an infinitely long carbon nanotube. 
The propagator between any two sites of an impurity-free nanotube is obtained and will be 
used to calculate the indirect exchange coupling between two magnetic adatoms. 

The tight-binding Hamiltonian for a carbon nanotube is diagonal when written in the 
basis |k, ±) defined as 



by a pair of indices (£, *) , where £ is the two-atom cell index and * = • or * = o represent 
the two inequivalent atomic sites of the hexagonal lattice. The quantity L represents the 



to the basis formed by two inequivalent carbon sites in the graphene structure. The band 
structure becomes 1-dimensional when quantization conditions on the wave vectors k x and/or 
k y are introduced. For achiral tubes, namely armchair and zigzag, the wave vectors k x and 
k y are quantized, respectively. In this case, the non-quantized wave vector runs continuously 
along the first Brillouin zone, whereas the quantized counterpart is only allowed to assume 
discrete values. For chiral tubes in general, a linear combination of k x and k y imposes the 
quantization condition. The resulting one- dimensional energy dispersions for the armchair 





The basis {£,*) represents atomic orbitals centred at a site located at R\ that is labelled 
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and zigzag nanotubes are given by^ 

Ei(k v : 

and 

Ei(k x ) = 



1 + 4cos (^) c os(^|) + 4c os 2 (^|) 



1 + 4cos( — - — )cos(— ) +4cos 2 (— ) 



(6) 



(7) 



respectively, where M depends on the geometry of the tube and is directly related to its 
diameter. The integer j = 0, 1,2, ...,M labels the different bands. The nearest-neighbour 
electronic hopping is denoted by 7 and is hereafter considered to be our energy unit. Notice 
that the bands are symetrically distributed around E = 0. 

In what follows we calculate an arbitrary matrix element of the Green function associated 
with both armchair and zigzag nanotubes, whose band structures are given by Eqs.tjHJ) and 
(J7J). Green functions for chiral tubes are also possible to be obtained within the same 
formalism, but the corresponding unit cells can be very large. The general formula for the 
single-electron Green function is 



Gip) = E 
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to 
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E+(k) uj-E_(k) 

where to is the energy. Introducing Eq.(j5J) into Eq.fjHJ), we obtain an expression for the 
Green function between any two sites \j, *) and of a graphene sheet. For simplicity, 

we consider the propagation between two equivalent sites in different cells. In this case the 
matrix element of the Green function is given by 



(t,*\g(u))\e,*) 



2L 
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A. Armchair nanotubes 

The sum over k in Eq.fJHJ) can be transformed into an integral over k y . The wave vector 
k x remains discrete and its values are k x = j^^^- I n addition, by taking into account that 
E-(k) = —E + (k), the expression above is rewritten as 

Jk.(R'-R' e/ ) 



AM-k 



■TV /a 

E / dk. 



lo e 



tt/o, Up' — E^(k Xj ky 



(10) 



For simplicity, in Eq. (fTU|) we have dropped the bra and ket representing the initial and 
ending sites. For each value of k x we have a corresponding one-dimensional band structure. 
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The integral above can be evaluated by extending k y to the complex plane and changing the 
integration contour from a straight line on the real axis to the boundaries of a semi-infinite 
rectangle in the upper half-plane whose base lies on the real axis between —ir/a and n/a. 
By determining the poles and their respective residues we can evaluate the integral for each 
band index k x . The poles, labelled qf, are solutions of the following equation 



t ± a \ 1 r r 3 n \ _l M ■ z/ 3 TTx-i ( ,,^ 
cos of- = — jcos ± \ — — sin \ , 11 

and the corresponding residues are 

R es r ±i = u e I (12) 

WjJ 47rMt 2 {cos(^f)sin(gff) + cos(of f) sin(g±f)}' V ; 

where Kf = ( J^ a )4 + (of)^y Eq. lfTTjl does not define the poles qf uniquely. By defining 
the cosine of pole qf but not its sine, the residue in Eq.(fT2|) is free to assume two distinct 
values. The residue depends on the correct sign of sin(g^|), which is obtained by imposing 
that the imaginary part of diagonal elements of retarded Green functions must be negative. 
The expression for Gt,t'{u>) then becomes 

ii? ± .(i?*-i?* ) 

yi,n ) ^ Mt2{cos( ^ f)sin(g ± f) + cos(g ± f)sin(g ± f)} - i ) 

This expression gives the matrix elements for the Green function of an armchair nanotube. 
Although the derivation of Eq. (jl3j) assumed equivalent sites in distinct cells, the expression 
is almost the same for inequivalent sites, except for a multiplying phase factor e~ % ^ K ). 

It is interesting to look at the physical significance of those poles. Each value of k x yields 
a corresponding k y and since they are obtained from the band structure for a given energy ou, 
they are just the coordinates of constant-energy surface along the axial direction. In other 
words, they indicate the wave vector qf with which electrons of energy to propagate along 
the axial direction of the nanotube. Many discrete values of k x have no real k y components, 
and therefore qf assumes imaginary values leading to evanescent contributions to the Green 
functions. 
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B. Zigzag nanotubes 



In the case of zigzag tubes the fcy-components are quantized whereas the k x run contin- 
uously over the Brillouin zone. The expression for the Green function becomes 

g f u ) = ^lly / dk x \ , 9/ -] , (14) 

e/K ! 2Mtt ^ J-^/aVs 1 uj 2 -El(k x ,k y y' V 1 

Analogously to the armchair case, the integral above can be evaluated by indentifying 
the poles and summing over the respective residues. The poles, labelled qj, are obtained 
through the following equation 

< "A ^-l-4cos 2 (ff) 

coster— ) = 2 -fj- (15) 

2 4cos( I f) 

and the corresponding residues are 

Res oJ = : j^, (16) 

47rMt 2 cos(f )sm( qj ^) 

where Kj = qj e x + (|p )e y . Finally, the expression for QiyiuS) is given by 

iuj e 1 **? 

GeAU) = ^ Mt2{cos( if ) sin(gff) + cos(qf § ) sin(gff)} ' 

This expression gives the matrix elements for the Green function of a zigzag nanotube. 

Despite having presented closed-form analytical expressions for the Green functions of 
armchair and zigzag nanotubes separately, common features in the expressions must be 
highlighted. The energy dependence is not only given by the energy on the numerator of 
Eqs.(JT3J) and (JHJ), but it is implicitly contained in the poles, as shown in Eqs. ffTTj) and (fT3j) . 
Likewise, the M-dependence is also implicit in the poles. Each term of the summation cor- 
responds to a band contribution to the single particle Green function. One major advantage 
of an expression for the Green function that does not depend on numerical evaluation, is 
that it allows a transparent analysis of the relationship between the electronic structure 
and the relevant parameters involved. For diagonal elements of the Green functions, j = f 
vanishes the argument of the exponential in the numerator of Eqs. (jl3J) and ()17j) and the 
imaginary part gives a direct expression for the density of states for both armchair and 
zigzag nanotubes. 
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IV. RESULTS AND DISCUSSION 



The indirect magnetic coupling can now be calculated by inserting the required matrix 
elements of Q into Eq.(@J) and evaluating the corresponding energy integral of Eq.(j3J). By 
changing the relative positions between the adatoms we can see how the coupling depends on 
their separation. Furhermore, by changing the nature of the adatoms we can also investigate 
how the coupling responds at the presence of different magnetic entities. We can also 
investigate how the coupling behaves for different types and sizes of nanotubes and finally, 
how it changes with temperature. 

A. Zigzag (metallic) 

We start by considering a metallic zigzag nanotube defined by the indices (6,0). The solid 
line of figure 121 shows the coupling AS(tt) as a function of the distance D between adatoms 
along the axial direction of the tube. The parameters were chosen to reproduce the adatom 
occupation of Co atoms. The first noticeable feature in this result is the monotonicity of 
the coupling. As a general physical principle, it is well know that the indirect coupling 
between magnetic impurities embedded in a metallic host oscillates as the impurities are 
moved further apart.— At first glance, one may think that this principle is not applicable to 
carbon nanotubes. However, as we show below, this is not the case, the lack of oscillations 
being the result of a commensurability effect. 

On a three-dimensional system, in an analogy to the de Haas van Alphen effect in metals, 
the oscillatory behaviour of the coupling is determined by the topology of the bulk Fermi sur- 
face (FS) of the host material.— More precisely, each FS extremum along a specific direction 
contributes with one oscillatory component of a definite period. This period is determined 
by the caliper associated with that extremum. For a CN, the Fermi "surface" is actually 
a discrete set of points. Inspecting the zigzag's Fermi points one sees that there is only 
one extremum, thus the coupling should oscillate with one single period. Nevertheless, the 
solid line in figure El is always negative indicating a preferential antiferromagnetic alignment 
between the adatom moments. For infinitely large separations, the coupling should vanish 
but it does so without sign changes or any oscillations. This is because the period coming 
from the caliper associated with that single extremum is actually commensurate with the 
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FIG. 2: Exchange coupling between Co adatoms as a function of the adatom separation on a zigzag 
(6,0) nanotube (solid line). The dashed line is the coupling when the nanotube is slightly doped, 
corresponding here to a small shift in its Fermi energy. The inset shows schematicaly the Fermi 
wave vector for pure (solid arrow) and doped (dashed arrow) tubes. All coupling values are given 
in units of the nearest-neighbour electronic hopping 7. 

nearest-neighbour distance in the axial direction. In other words, the expected oscillations 
are hidden by the fact that the coupling is being probed at intervals that are of the same 
lenght as the unit cell separation. To prove this point, we introduce a slight shift in the 
Fermi level that would correspond to either doping the nanotube or to the action of a gate 
voltage affecting the occupation of the system. In this case, the Fermi surface of the nan- 
otube no longer consists of isolated points at the corners of the Brillouin zone. The points 
in reciprocal space are now located at the intersection between small circular-like curves 
centred at the hexagonal corners of the Brillouin zone and the discrete quantization lines 
along the k y direction, as shown in the inset of figure EJ The inset displays the wave vectors 
that contribute to the oscillatory coupling, which clearly shows that they no longer coincide 
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with the dimensions of the Brillouin zone. A small shift in the Fermi level (AEp = O.O67 in 
this case) is sufficient to destroy the commensurability effect in the coupling, as shown by 
the long-period oscillatory curve (dashed line) of figure 121 

The lack of oscillations in the indirect coupling is sufficiently robust to appear with other 
magnetic adatoms. In fact, by changing the parameters such that the d-band occupation of 
the adatoms corresponds to that of Ni atoms, we find a similar behaviour for the coupling, 
the only difference being the overall sign. This is illustrated by figure El in which the full line 
depicts the calculated coupling between two Ni adatoms as a function of their longitudinal 
separation. Rather than an antiferromagnetic alignment between the moments, in this case 
we find that they tend to align in a ferromagnetic fashion. Likewise, the monotonicity of 
the coupling remains unnaffected by varying the tube diameter. This is confirmed by the 
dashed and dot-dashed lines in figure H3 in which the coupling between Ni adatoms is plotted 
as a function of their separation for a (9,0) and a (12,0) tube, respectively. The inset also 
shows that the coupling decays with the tube diameter as 1/R, where R is the nanotube 
radius. This is hardly surprising, bearing in mind that the density of electronic states of a 
nanotube has a similar behaviour and that the coupling is the result of an indirect coupling 
mediated by the nanotube electrons. 

B. Zigzag (semiconductor) 

For the case of semiconducting tubes, the lack of extended states at the Fermi level 
supresses the long-range character of the indirect coupling, even in 1-dimension. In this case 
the coupling decays exponentially with the adatom separation, as shown by the full line 
in figure |U for a (31,0) nanotube. The choice of such a large diameter tube is justified to 
illustrate the fact that a minute change in the occupation is capable of moving the Fermi 
level into the zone of extended states and producing a remarkable effect in the coupling. 
The dashed line shows the coupling for a slightly doped nanotube (AEp = O.O67) displaying 
similar behaviour to the metallic case, where a long-period long-ranged oscillatory coupling is 
observed. From a virtually uncoupled configuration, the coupling between adatom moments 
can be switched on by changing the occupation of the nanotube. Consequently, by using 
gate voltages to control the Fermi level of the system, in principle one can reversibly turn 
magnetically uncoupled adatoms into fully aligned magnetic moments, a promising feature 
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FIG. 3: Exchange coupling as a function of longitudinal separation D between Ni adatoms on (6,0) 
(solid line), (9,0) (dashed line) and (12,0) (dot-dashed line) nanotubes. The inset shows how the 
coupling amplitude decays as a function of tube diameter. All coupling values are given in units 
of the nearest-neighbour electronic hopping 7. 



for use in spintronic devices. 



C. Armchair 



For the armchair tubes, there are two contributory wave vectors extracted from their 
Fermi surface, and as such there should be two oscillatory components in the exchange 
coupling. These are identified in the inset of figure together with the coupling as a 
function of the adatom separation in the main graph. The full line corresponds to the 
coupling evaluated at T = and the dashed line shows the coupling at room temperature. 
Once again, one of the periods is commensurate with the interatomic distance, which makes 
the coupling to oscillate with a non-commensurate period superimposed to a steady decay. 
The coupling at finite temperature tends to decay faster, which can be explained by the 
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FIG. 4: Exponentially decaying exchange coupling between Co adatoms on a semiconducting (31,0) 
nanotube (solid line). A small amount of doping is enough to populate the tube's conduction band 
and give origin to an oscillating, long-ranged exchange coupling (dashed line). All coupling values 
are given in units of the nearest-neighbour electronic hopping 7. 

fact that the thermal energy partially suppresses the propagators between adatoms, thus 
reducing the efficiency with which the electronic carriers transport the magnetic information 
back and forth. 

Concerning the rate of decay of the coupling as a function of the adatom separation, 
all examples above decay slowly in the asymptotic limit (D >> 1). It is known that in 
one-dimensional structures the coupling decays approximately as 1/D for D » 1. This is 
what we observe approximately when impurities embedded in the nanotube are considered, 
instead of adatoms. As the tube diameter increases, the system approaches a 2-dimensional- 
like structure and the coupling tends to decay a little faster. This is displayed in Figure 
El showing the approximate decaying rate as a function of the tube diameter for armchair 
structures. The coupling between adatoms has a more complicated asymptotic behaviour. 
For most of the range of D values we investigated it was not possible to identify a unique 
decay rate. In the regions where the decay rate is well-defined it may be even slower than 
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FIG. 5: a) Exchange coupling as a function of longitudinal separation D between two Co adatoms 

on a (6,6) armchair nanotube. The inset shows the wave vectors contributing to the coupling with 

different oscillation periods, b) Amplitude of the exchange coupling as a function of temperature. 

All coupling values are given in units of the nearest-neighbour electronic hopping 7. 
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FIG. 6: Exchange coupling A£(ir) multiplied by the distance between impurities D, as a function 
of D, for a) (6,6), b) (24,24) and c) (48,48) armchair nanotubes. It is evident that in the (6,6) and 
(24,24) nanotubes the coupling decays as l/D, whereas in the (48,48) tube it decays faster. All 
coupling values are given in units of the nearest-neighbour electronic hopping 7. 

l/D. Differences between the coupling decay rates for adatoms and embedded impurities 
have been reported in other systems^. 

It is worth mentioning that such a slow decaying rate in the coupling leads to a long- 
ranged magnetic interaction that, when combined with the lack of oscillations introduced by 
the commensurability effect, can be potentially used to induce ordered magnetic structures 
on the surface of a carbon nanotube. 



V. CONCLUSIONS 

In summary, we have presented expressions for the indirect exchange coupling between 
magnetic adatoms mediated by the conduction electrons of carbon nanotubes. It is obvious 
from the expressions that the coupling, written in terms of a few matrix elements of the 
single-particle Green functions, depends on the electronic carriers transporting the magnetic 
information back and forth to both adatoms. We show that, contrary to expectations, 
we find a monotonic behaviour in the coupling between magnetic adatoms adhered to the 
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surface of metallic zigzag nanotubes. Rather than an intrinsic property, we have shown 
that the lack of usual oscillations results from a commensurability effect. In other words, 
the oscillations are hidden by the fact that the contributory wave vectors have the same 
size as the Brillouin zone, forcing the periods to coincide with the atomic spacing. In this 
commensurate regime the coupling does not change sign but decays monotonically. The 
overall sign of this monotonic coupling depends on the nature of the adatoms attached 
to the nanotube. In the case of positive coupling, a FM alignment between the moments 
persists as we move the adatoms apart, where the coupling magnitude decays rather slowly as 
approximately 1/D. Such a long range interaction may be responsible for inducing magnetic 
order in the nanotube. 

The coupling across semiconducting tubes is exponentially damped, which results from 
the lack of extended states to transport the magnetic information between the adatoms. 
By moving the occupation into the conduction band we show a dramatic change in the 
coupling, departing from a configuration with random alignment between magnetic moments 
into another that is highly ordered due to the predominantly positive long-ranged long- 
period oscillatory coupling. We argue that a controlled gate voltage can potentially tune 
the ocuppation in such a way that the magnetic interaction can be switched on and off in a 
reversible fashion. 

Oscillatory coupling is found on armchair metallic tubes. In this case, two periods are 
identified, one commensurate and another incommensurate with the lattice. The result- 
ing coupling consists of a monotonically decaying contribution superimposed to a 1/D- 
modulated oscillatory coupling. Concerning the dependence of the coupling on other pa- 
rameters, its diameter dependence follows the l/i?-behaviour displayed by the density of 
states. As expected, the coupling decays with temperature but not as fast as the coupling 
between moments embedded in 3-dimensional structures. 

Finally, it is worth emphasizing that by controlling the alignment of magnetic structures 
attached to the surface of carbon nanotubes one can potentially select the way in which 
the systems will respond to magnetic excitations. Moreover, bearing in mind that the 
transport properties of magnetically doped structures are sensitive to the alignment of their 
magnetizations, the understanding of how separate moments are coupled across nanotubes 
could be a major step forward in the search for and control of spin-valve effects in those 
structures. 
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